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The problem of estimating an unknown force driving a linear oscillator is revisited. When using linear mea- 
surement, feedback is often cited as a mechanism to enhance bandwidth or sensitivity. We show that as long as 
the oscillator dynamics are known, there exists a real-time estimation strategy that reproduces the same mea- 
surement record as any arbitrary feedback protocol. Consequently some form of nonlinearity is required to gain 
any advantage beyond estimation alone. This result holds true in both quantum and classical systems, with non- 
stationary forces and feedback, and in the general case of non-Gaussian and correlated noise. Recently, feedback 
enhanced incoherent force sensing has been demonstrated [Nat. Nano. 7, 509 (2012)], with the enhancement 
attributed to a feedback induced modification of the mechanical susceptibility. As a proof-of-principle we ex- 
perimentally reproduce this result through straightforward filtering. 



Micro and nano-mechanical oscillators are capable of ultra- 
sensitive force measurement, allowing precision spin, charge, 
acceleration, and field sensing G]^. It is well known that 
linear feedback control can improve the performance of non- 
linear mechanical sensors |'5^,'6|. For example, in non-contact 
atomic force microscopy linear feedback is commonly used 
to stabilise the tip-surface separation, thereby avoiding col- 
lisions and suppressing frequency drifts due to short range 
forces such as van der Waals forces QUI- Since feedback 
control modifies the response of an oscillator to environmental 
forces it also appears attractive as a technique to enhance the 
performance of linear sensors. For instance feedback cooling 
allows the suppression of thermal noise |I9]|T3, while feed- 
back tuning of the spring constant can provide increased me- 
chanical response at the signal frequency 1 11 1. However, such 
precision enhancement is prohibited for linear processes with 
stationary linear feedback and uncorrected Gaussian noise by 
the well-known principle of neutrality in linear control theory, 
which states that the accuracy with which the oscillator posi- 
tion can be determined is independent of feedback lfT2l[T3l . 

Non-stationary processes, non-Gaussian noise and non- 
linear estimation strategies are each found in a range of lin- 
ear oscillator-based force sensors. Stroboscopic measurement 
of impulse forces |fT4l , and variance estimation of incoher- 
ent forces as applied in bolometry lITSI ITSll . are two relevant 
examples. Linear feedback cooling has been proposed as a 
means to enhance precision in both cases lfT4l[T7l . and exper- 
imentally demonstrated in the latter ifTTl . However, neither 
proposal identifies an optimal estimation strategy. This leaves 
unresolved the important question of whether the same, or im- 
proved, sensitivity enhancement might be achieved without 
feedback by applying a better estimation strategy. 

Here we present a straightforward theoretical approach 
which shows that, even in the presence of non-Gaussian or 
coiTelated noise and non-stationary processes, a real-time es- 
timation strategy always exists that reproduces the same mea- 
surement record as any arbitrary linear feedback protocol (see 
Fig. The theory applies to both quantum and classical 
oscillators, and to the intrinsically non-linear problem of vari- 
ance estimation IfTTl . It ultimately provides a clear set of min- 
imum requirements for feedback to provide any advantage to 



force sensing over that possible with estimation alone. Es- 
sentially, some form of nonlinearity is required in the physi- 
cal system, which may arise from the measurement process, 
feedback loop, signal, or from the oscillator itself. When no 
non-linearities are exhibited, the theory yields a filter which 
allows the estimate that would be obtained with feedback to 
be determined causally from the measurement record without 
feedback. This precludes the possibility of any additional sen- 
sitivity enhancement from feedback in either of the examples 
discussed above lll4HT7ll . or indeed feedback improved band- 
width in linear force sensors ITSI . 




FIG. 1 : (Color online). A Conceptual diagram comparing optome- 
chanical force sensing with feedback and filtering B Experimen- 
tal schematic. Red (dark grey): fiber interferometer; solid green 
(light gray): electrical components for feedback stabilization; dashed 
green: electrical signal applied to the microtoroid for application of 
the incoherent force and feedback. FPC: fiber polarization controller. 
BPF: bandpass filter. MZM: Mach-Zender Modulator. PM: phase 
modulator. PI: proportional-integral controller 

To validate the theoretical results, we experimentally repro- 
duce the effect of feedback enhanced sensitivity achieved in 
Ref. IfTTl . but replacing feedback with causal filtering. This 
demonstrates that the sensitivity enhancement achieved in that 
experiment was not due to a feedback-induced change in cou- 
pling of the oscillator to its the environment, as is suggested 
in Ref. [1T|. Rather it arises from an intrinsic bias in the esti- 
mator used towards low quality, and thus feedback cooled, os- 
cillators. By clearly demarcating the circumstances in which 
feedback may advantage force sensing over estimation alone. 
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our results both clarify a significant ambiguity in the op- 
tomechanics and force sensing communities, and contribute 
towards simplifying experimental implementations of ultra- 
precise force sensing with linear oscillators. 

The evolution of a mechanical oscillator can be described, 
in both classical and quantum regimes l,19 j by the equation of 
motion 



feedback force is 



m\x + yx + Q,^,.x] ^F^{t, x) + F.^i{t, x) 



(1) 



where m, y, and Q„, are the mass, damping rate, and reso- 
nance frequency, respectively. For compactness the combined 
force Fj^{t,x) - Fj{t,x) + F^{t,x) + F\,J,t,x) is used, where 
F-iit, x) is a thermomechanical force due to the coupling of 
the oscillator to its environment, F^it, x) is the signal force 
and Fba(f, x) is a backaction force due to the act of measure- 
ment. F^ait, X) is an actuation force used for feedback where 
x(t) - x(t)+N(t, x) is the instantaneous measurement record of 
the oscillator position x, and A^(f, x) is the measurement noise 
which maybe correlated to the backaction noise. In general, 
the forces and measurement noise can all have non-stationary 
dynamics, non-Gaussian noise and non-linear dependence on 
the oscillator position. They can each be linearized by Taylor 
expanding around the mean position of the oscillator x and re- 
taining only zeroth and first order terms (See Supplementary 
Information). The zeroth order terms are independent of fluc- 
tuations in the oscillator position, forcing, and measurement 
noise; and only serve to deterministically shift the mean po- 
sition of the oscillator. The first order terms are each linearly 
dependent on only one source of fluctuation; and either act to 
modify the mechanical susceptibility or introduce incoherent 
driving. Higher order terms introduce nonlinearities and insta- 
bilities which can give rise to detrimental effects such as sat- 
uration and nonlinear dissipation [5^, 20]. By neglecting these 
higher order terms, we restrict the analysis to the most general 
linear oscillator experiencing linear feedback. It is important 
to note that the deterministic shift in mean oscillator position 
due to the zeroth order terms can affect the force sensitivity; 
for example, by shifting a cavity optomechanical system onto 
optical resonance. However, since this is deterministic, and 
known a priori, an equivalent displacement may be made to 
the oscillator without feedback by applying a known external 
force, as depicted in Fig.[T|\ (see Supplementary Information 
for details). 

To simplify the analysis and present results most relevant 
to our experiments, in the main text of the Letter we consider 
the common scenario where the mechanical oscillator's sus- 
ceptibility is only modified by the feedback force, and there- 
fore drop the first order susceptibility modifying terms in the 
other forces. We further assume that all processes involved 
are stationary. These specific assumptions are not necessary 
for our conclusions, which hold for the most general linear 
case, including non-stationary processes and first order terms 
(see Supplementary Information). Under these assumptions, 
the combined force is given by F^nit, x) = F^nit, x) and the 



F^ci(t, x) 
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git - T)x{T)dT 
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where x(t) - x{t) - N(t, x) and g(t - r) is the stationary feed- 
back kernel describing the filter applied to the measurement 
record. Enforcing causality, namely g(t - t) = for t > t, 
simplifies Eq. (j2|i into the convolution Fact(f, x) - g{t) * x(t). 
Substitution into Eq. ([T]i and Fourier transforming then gives 



x(Q) - x(^) [^m(fi, x) + g(Q)i(Q)] 



(3) 



where xi^T^ = '«J^^m - ^'^ + ij^] 

is the intrinsic mechani- 
cal susceptibility 1281 . The oscillator position without feed- 
back can be simply obtained by omitting the feedback force 
g(Q.)x(Q.) from Eq. ([3]l, xo(Q) = ;t'(Q)F„(Q, xq), where the 
subscript is used to distinguish from the feedback case. 

As discussed earlier, an external force may be applied to 
equate the mean positions of the oscillator with and without 
feedback. In the case of cavity optomechanics this amounts 
to ensuring identical cavity detunings when experiments are 
initiated. With x = xq, the common forcing terms with and 
without feedback F^(Q., x) and F^(n, xq) are identical. Sub- 
stituting for x(Q) and xo(fi) in terms of their respective mea- 
surement records (eg. x(Q.) - x(D.) - N(D., x)) then gives a 
completely deterministic equation relating the time domain 
measurement records that would be achieved with and without 
feedback 
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where h(Q.) - x' Ix is the ratio of the modified mechanical sus- 
ceptibility to the intrinsic mechanical susceptibility. There- 
fore the exact position record that would be obtained using 
stationary linear feedback can be retrieved straightfowardly 
by applying the filter h{D) to the position record without feed- 
back. This precludes enhancements of both sensitivity ifTTll 
and bandwidth ifTSl beyond that achievable with estimation 
alone. Since no constraints are placed on the statistics of the 
driving forces or measurement noise, this result is valid even 
for non-Gaussian noise or if correlations exist between mea- 
surement and process noise, such as those induced by quan- 
tum backaction. Furthermore, since it applies directly to the 
measurement records, rather than a specific parameter esti- 
mation process based on them, it holds for both linear and 
non-linear estimation processes. We show in the supplemen- 
tary information it can be generalised to include linear non- 
stationary forcing and feedback as well as modifications to 
the mechanical susceptibility due to effects such as optome- 
chanical dynamical backaction IfTOl |2T1 l22l . In this case, 
the required filter is more complex and is, in general, non- 
stationary, but remains causal. Our results are valid in both 
the quantum and classical regime and rigorously prove that no 
force sensing advantage is provided by linear feedback onto a 
linear oscillator with known dynamics. Consequently, nonlin- 
earities are a minimum prerequisite for feedback to improve 
force sensing beyond estimation alone. 
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Recently, enhanced incoherent force sensing was experi- 
mentally demonstrated 111 71 by stationary feedback cooling 
of a linear oscillator. However, as shown here, no sensitiv- 
ity enhancement is obtained from this method over estimation 
alone. The exact filter equivalent to the feedback cooling used 
in Ref. IfTTl is obtained by substituting g{Q.) = -/myOgf into 
Eq. (j4]), where gf represents the filter's unitless gain. This fil- 
ter, denoted hd^), effects the causal map xq i-> x. We demon- 
strate this experimentally here in a similar system to that of 
Ref 1 17 1 consisting of a microtoroidal cavity optomechani- 
cal system. An intrinsic mechanical mode of the microtoroid 
is used to transduce an incoherent electrostatic gradient force 
applied by a nearby electrode |23|. A whispering gallery opti- 
cal mode of the microtoroid is used to read out the mechanical 
motion and thereby determine the variance of the incoherent 
force. 

Our experimental setup is shown in Fig. [TJi. A shot-noise 
limited fiber laser at 1550 nm was evanescently coupled into 
the whispering gallery mode of the microtoroid using a ta- 
pered optical fiber. The microtoroid had major and minor di- 
ameters of 60 /vm and 6 fim respectively with a 26 fim under- 
cut. The mechanical motion of the microtoroid, which induces 
phase fluctuations on the transmitted light, was detected inter- 
ferometrically by beating with a bright 3.5 mW optical phase 
reference followed by shot-noise limited homodyne detection. 
We actively stabilize the toroid-taper separation using an am- 
plitude modulation technique [!24j that maintains a constant 
coupling rate into the optical cavity. Pound-Drever-Hall lock- 
ing was used to lock the laser frequency to the optical reso- 
nance, which had an intrinsic quality factor of Qo =2.6 x 10^. 
A 50/50 tap-off' after the microtoroid was used to detect the 
cavity transformed amplitude and phase modulation, provid- 
ing the error signal for the optical frequency and taper-toroid 
separation locks. The interferometer was locked midfringe 
via a piezo actuated fiber stretcher that precisely controls the 
optical path length in one arm. 

The measurement record is acquired from the homo- 
dyne signal by electronic lock-in detection where demodu- 
lation of the photocurrent at the mechanical resonance fre- 
quency allows real time measurement of the slowly evolv- 
ing quadratures of motion, denoted /(f) and Q{t) where x{t) - 
/(f)cos(Q,„r) + Q{t)sm(Q.,„t). Fourier analysis reveals a me- 
chanical power spectra with peaks corresponding to micro- 
toroid mechanical resonances. Fig.[2|\ (red) shows the room 
temperature Brownian motion of a mechanical mode with a 
signal-to-noise ratio of 37 dB and a fundamental frequency, 
damping rate and effective mass of Qm = 40.33 MHz, y - 
23 kHz and m^ff = 0.6 fig respectively. The absolute mechan- 
ical displacement amplitude was calibrated via the optical re- 
sponse to a known reference phase modulation [25 j . 

As shown earlier, applying the filter h^i^) to the mea- 
surement record without feedback should retrieve an identi- 
cal measurement record to that obtained with feedback. Ap- 
plying this filter to the measurement record it is possible to 
mimic feedback cooling as shown in Fig. |2j\ (green) where 
the filter gain is varied from 2 to 150. Extending the gain 



1 10-^^ 



10-' 





I a 0.5 

A I ° 2 
A 


A 




/ \ 20 






Jr ^ -20 








-1 0,1. 
5m/2it (10 Hz) 




0.2 0.5 1 2 

Averaging time (ms) 



10" 10' 
Filter Gain 



FIG. 2: (Color online). A Displacement spectrum at room tempera- 
ture (red) and with filtering (green) at gain 2,8,34,72 and 150. Inset 
top: displacement spectrum with feedback cooling (blue) and filter- 
ing (green) with gain I. Inset bottom: percentage difference between 
filtered and cooled spectrum with equivalent gain. B Force sensitivity 
as a function of averaging duration for thermal (red points), feedback 
cooled (blue points) and filtered spectra (green points) at gain 1,2.4,5 
and 10. Red line: fit to inverse quartic dependence of the force sen- 
sitivity without filtering or feedback. C Force sensitivity versus filter 
gain after 1 ms of averaging (green circles) showing good agreement 
to theory (dashed line). 



beyond > 20 the mechanical spectrum inverts and exhibits 
squashing, a well known characteristic of high gain feedback 
cooling |23|. To confirm the equivalence of the measurement 
record obtained via feedback and filtering we implement feed- 
back cooling by applying the homodyne photocurrent to the 
toroid through an electrode which generates strong electrical 
actuation of the mechanical motion through electrical gradient 
forces |26|. This allows the mechanical mode to be cooled 
from room temperature by a factor of 2. The upper inset in 
Fig|2j\ shows feedback cooling (blue) and equivalent gain fil- 
tering (shaded green); with the fractional difference between 
the feedback and filtering spectra showing no statisticafly sig- 
nificant difference (lower inset). 

An estimate of the variance of an incoherent force applied 
to an oscillator may be obtained by determining the oscilla- 
tor's energy 1 171. After averaging time, r, the estimate of the 
energy is given by = 1 /t dtl(t)^ + Q{t)^. To calculate the 
ensemble average {Ej) and the standard deviation cr^{f) of 
the energy estimate multiple independent measurements are 
made for each r. Following Ref. [17] the energy estimate can 
be translated into an estimate of the magnitude of the force 
with a sensitivity given by 6F{t)^ - o-Eir)/ J^^ dD.\x' ■ It 
is important to note that this estimation process is not neces- 
sarily optimal. In the case where r > 1/y the force sensitivity 
scales as (yr)"'^"* which appears to motivate the use of feed- 
back cooling to increase the mechanical decay rate y IJ7i . 
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Figure (red points) shows the inverse power-law depen- 
dence of the force sensitivity on averaging duration, t, for our 
experiments with only thermal driving. As predicted by our 
theory, by applying the filter h^iO.) to the thermal data it is pos- 
sible to enhance the force sensitivity in the same way as feed- 
back cooling. This is shown in Fig.|2jj (green) where increas- 
ing the filter gain, gf, provides a clear improvement in sensi- 
tivity while consistently maintaining the predicted power-law 
dependence with respect to averaging time. Figure [2p (cir- 
cles) shows the force sensitivity as a function of filter gain 
taken for a fixed averaging duration of t = 1 ms. For gains 
below = 20 the measured sensitivity agrees with the theo- 
retical fit. At higher gains it is degraded due to squashing of 
the mechanical power spectrum which suppresses the signal 
and allows shot noise to dominate. 
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FIG. 3: (Color online). A Normalised energy estimate at room tem- 
perature (red), and with addition of incoherent driving (green). Solid 
lines: mean; dashed lines: one standard deviation error bounds. B 
Averaging time required to resolve incoherent signal versus filter 
gain (green circles). Dashed line: theory 



To demonstrate the improvement in sensitivity achievable 
via filtering we apply a small incoherent electrostatic gradient 
force to the microtoroid with a magnitude of approximately 
7% of the thermal energy. The ability to resolve this force 
against the thermal noise depends on the averaging time. Only 
when the standard deviation of the energy estimate is smaller 
than the strength of the signal can the incoherent force be re- 
solved. The convergence of the thermal energy estimate with 
increasing averaging time is shown in Fig.[3]\ (red). With the 
addition of the incoherent signal the ensemble average is in- 
creased without affecting the error bounds as shown in Fig.|3]3 
(green). At 3ms the error becomes comparable to the energy 
separation and the applied incoherent force is resolved. If the 
filter hc{Q.) is applied the force sensitivity is improved and the 
time taken to resolve the applied force decreases as shown in 
Fig.[3ji. The required averaging time decreases as the estima- 
tion gain is increased in good agreement with theory. At gains 
gf > 20 the averaging time increases again due to inversion 
of the mechanical spectrum and suppression of the signal rel- 
ative to shot noise. The inflection point in Fig.|3ji and Fig.|2p 
shows that even though thermal noise dominates shot noise by 
orders of magnitude at the peak of the mechanical susceptibil- 
ity, in incoherent force sensing it is shot noise that determines 
the ultimate sensitivity limit. 



The experiments presented here show that feedback and fil- 
tering allow equivalent enhancement in incoherent force sen- 
sitivity. In this context, the results of Ref. ITTll . can be nat- 
urally understood as a consequence of an intrinsic bias that 
arises when using the oscillator energy to perform incoherent 
force estimation. Near-resonant spectral components of the 
incoherent force drive the oscillator more strongly, and are 
therefore over-represented in the measurement. As a result, 
even though feedback only applies a reversible transforma- 
tion to the measurement record, the force sensitivity appears 
to improve with increasing oscillator linewidth as the estima- 
tion becomes more balanced. The stroboscopic feedback en- 
hanced force sensing scenario proposed in [14] is similarly 
biased. In that case, measurements prior to application of the 
signal force are used to pre-cool the oscillator. This improves 
its initial localization in phase space and, thereby, the capac- 
ity to resolve displacements due to external forces. However, 
the existence of this prior measurement record is not taken 
into account when calculating the sensitivity of measurements 
without feedback. Filtering it appropriately allows equivalent 
localisation to feedback cooling, though offset from the origin, 
and achieves identical sensitivity. These examples illustrate 
the main result of this Letter that, for a linear oscillator, any 
sensitivity enhancement arises not through the action of feed- 
back, but rather through measurement and estimation alone. 

In summary, we have theoretically shown that for linear 
oscillators neither stationary nor non-stationary linear feed- 
back provide any force sensing enhancement over that possi- 
ble with estimation alone. This holds true in both the quan- 
tum and classical regime and in the presence of non-Gaussian 
noise. As a demonstration we have shown experimentally that 
detection of a stationary incoherent force can be enhanced 
equally effectively via estimation as with feedback cooling. 
These results may have broad relevance to many scientific and 
engineering communities, particularly those associated with 
mechanical sensors limited by thermomechanical noise. Ac- 
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Supplementary Information 



expanded without loss of generality, as 



Here we show that the theoretical result presented in the 
main text of the Letter can be generalised to include linear 
non-stationary forcing and feedback as well as modifications 
to the mechanical susceptibility due to effects such as op- 
tomechanical dynamical backaction IfTOl 1211 l22l . Due to non- 
stationarity, the forces cannot be expressed as a simple convo- 
lution and hence requke a more complex theoretical approach. 

The generalized optomechanical system considered here is 
shown schematically in Fig.[T[\. The evolution of the mechan- 
ical oscillator can be described, in both classical and quantum 
regimes |fT9l by the following equation of motion 



i[x + yx + nix] ^FUt, X) + F^^tit, x) 



(5) 



where m, y, and Q,„ are the mass, damping rate, and resonance 
frequency of the mechanical oscillator, respectively. For com- 
pactness the combined force - Fj(t, x)+Ff^{t, x)+F^,.^{t, x) 
is used, where Fjit, x) is a thermomechanical force due to the 
coupling of the oscillator to its environment, Fs(f, x) is a gen- 
eral signal force and F^^it, x) is a backaction force due to the 
act of measurement. F^^tit, x) is a linear actuation force where 
x{t) = x{t) + N(t, x) is the instantaneous measurement record 
of the oscillator position x, and N(t, x) the measurement noise 
which maybe correlated to the backaction noise. In general, 
the forces and measurement noise can all have non-stationary 
dynamics, non-Gaussian noise and non-linear dependence on 
the oscillator position. 



x+N 




FIG. 4: (Color online). A General theoretical schematic outlining 
the emergence of forces and feedback. B Experimental schematic. 
Red (dark grey) indicates the fiber based interferometer; Solid green 
(light gray) indicate the electrical components for feedback stabiliza- 
tion and dashed green represent the electrical signal applied to the 
microtoroid for feedback and application of incoherent force. FPC: 
Fiber polarization controller. BPF: Bandpass filter. MZM: Mach- 
Zender Modulator. PM: Phase modulator. PI: Proportional-integral 
controller 



The combined thermal, signal, and backaction force may be 



Fm(t,x) = I dTgm(t,T,x) +n^(t,T,x)^i-n(T) (6) 

Jo 

= I drg^it, T, x) + g'^it, T, x)(x(t) - x) 
Jo 

+nUt,T,x)^UT) (7) 

= I dTGUt,T,x) + g'^{t,T,x)(x(T)-x) (8) 
Jo 

where the first term in Eq. (|6]l, gm(t, t, x), is a generalized gain 
kernel that allows an arbitrary coherent response at time f, de- 
pendent on the full prior history of the oscillator position x. 
The second term includes all incoherent forcing with ^m{T) 
being a unitless white-noise Wiener process and «m(r, t, x) a 
memory kernel that permits non-Markovian, non-stationary 
and position dependent noise. The actuation force can be ex- 
pressed similarly by replacing x with x in Eq. (|6]l. In general, 
both terms in Eq. (j6| can introduce nonlinearities and instabil- 
ities which can give rise to detrimental effects like saturation 
and nonlinear dissipation |5 20|. Here, we wish to show that 
no advantage is possible from feedback in the case of a hn- 
ear oscillator and linear feedback. Consequently, we Taylor 
expand both memory kernels gm(t, t, x) and t, x) about 
the mean position of the oscillator x and neglect nonlinear 
terms. This gives gm{t,T,x) = g^(t,T,x) + g'^{t,T,x)(x - x) 
where g'^{t,T,x) = t,x)/(?x|j.=^ for the coherent mem- 

ory kernel. Since the incoherent memory kernel is multi- 
plied by the noise term ^m(T"), nonlinearity occurs even in the 
second term of the expansion, so that for a linear oscillator 
nm(f, T, x) = n,n(f, T, x). Substituting these expressions into 
Eq. (j6| results in Eq. (|7]i, which can be further simplified into 
Eq. Jsjl by grouping all terms without dependence on the in- 
stantaneous displacement of the oscillator from its mean posi- 
tion into a single term Gm(f, t, x)-gm(t, t, x)+nm(t, t, x)^m(T) 
containing both noise and deterministic forcing. 

The actuation force Fact can be similarly Taylor expanded 
into its component force as 



F^ctit, x) 



-f 

Jo 



drG^ctit, T, x) + g'^^iit, T, x)(x(t) - x). (9) 



where x(t) = x(T)-N(t,x). Here, the term g^j,j(f,T,x)(x(T)-x) 
provides linear feedback modifying the susceptibility of the 
oscillator and driving its motion with white measurement 
noise N(t,x). Taking the case of ideal feedback where no 
noise is introduced by the feedback loop itself, Gact(f, x) = 
^act(f, x), acts only to shift the mean position of the oscilla- 
tor independent of the measurement record. Substituting the 
expressions for the actuation and combined forces into Eq. (jSj) 
and Fourier transforming gives 

x(Q) = ;r(^^)-^,-n| j| dTC^it, r, -x)+g'j,t, r, x)(x(t)-x) 

+ Gact(f, T, X) + g'^,{t, T, x)(x(t) - x)| (10) 
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where ^t^n represents the Fourier transform and x(^) ^ - 
m ^Q.^^ - + /yilj is the intrinsic mechanical susceptibility. 
It should be noted that the Fourier transform can only be taken 
if the oscillator is intrinsically stable. Consequently, scenarios 
such as regenerative amplification where feedback gives rise 
to exponential growth in displacement and eventually triggers 
a nonlinear response Il22l |27| are excluded from the analysis 
here. 

If feedback is not applied, both the linear feedback term 
and the feedback induced deterministic displacement are re- 
moved. In principle, the resulting change in mean oscillator 
displacement could modify the force sensitivity. For example, 
in cavity optomechanics, the displacement could move the op- 
tical cavity toward resonance, modifying the optical power 
level in the cavity and through this the sensitivity. However, 
since the feedback induced displacement is deterministic and 
known, an equivalent displacement may be achieved without 
feedback by applying an external force. In the case of cav- 
ity optomechanics, operationally, this corresponds to ensuring 
the cavity is detuned from resonance equivalently in both the 
feedback and non-feedback scenarios at the beginning of the 
measurement. Applying this external force so that x = x^, 
the oscillator position without feedback can then be simply 
obtained by omitting the feedback force g'^^^it, t, x){xiT) - x) 
from Eq. (fTOll 



xaiO.) = xi'^)^.->a{ dTGUt, t, x)+g'^(t, t, x){xq{t) - x) 



-l-Gact(r, T, X) 



(11) 



where the subscript is used to distinguish from the feedback 
case. 

Subtracting Eq. (10 1 from Eq. (Ill eliminates the com- 
mon terms G^^i{t,T,x) and Gn,(f, r, i:) and, after substituting 
x(Q) - x(D.) - N(Q., x), results in a completely deterministic 
equation relating the time domain measurement record with 
and without feedback 

m - j| dT[g:st\T,m8:jt\TMx(T)-x]l 

= io(f)-^oi, |;r(i^)^,'^n|j|«'r£(f',T,i)(io(T)-i)|l2) 

where the dummy variable f' has been introduced to distin- 
guish the Fourier transform from its inverse. Eq. ( 12 1 gives 
a non-trivial relationship between the measurement records 
with and without feedback in the presence of non-stationary 
forcing and feedback, non-Gaussian noise and correlations be- 
tween measurement and process noise. However, for feedback 
and filtering to be equivalent, a causal filter must exist that 
maps xo 1-^ X. To determine the existence of such a filter it is 
necessary solve for x(t) as a function of xi)(t) in Eq. (fT2]i. 

Without loss of generality we may simplify Eq. ( |12| l by 
choosing our position coordinate such that the mean displace- 



ment is zero x = 0. The second term on the LHS of Eq. ( 12 1 
may then be expanded as 



dt' dT[g'Jt',r,x)+gUt',T,x)\x(T)xit-t') (13) 

CO Jo 
fOO V /-*oc 

= I t/T I dt'[g'Jt',T,x)+g'^Jt',T,x)]x(t-t' 

Jo [J-oo 

- I dr {hm(t, T, x) + /iact(f , T, x)) x(t) 
Jo 



X(T)(14) 

(15) 



where the time shift property of Fourier transforms 
£^ dQ.e'^'-'''">X(n) = x(t-t') has been used to obtain Eq. ([ojl. 
Eq. (14 1 is found by enforcing causality, namely g'jit', t) - 
when T > t', allowing the upper bound of the integral over 
T to be extended to infinity and the order of integration to 
be rearranged. The term contained in the large square brack- 
ets of Eq. ( [T4| is a combined transfer function for the system 
which can be consolidated into terms denoted h^(t, t, x) and 



/!a(.t(r, T, x) that are zero when t > f, resulting in Eq. ( 15 1. Ap 



plying the same procedure to the RHS of Eq. ( 12i simplifies 
the equation into a Fredholm equation of the second kind 



x(t)- f 
Jo 



dT (hin{t, T, x) + h^ctit, T, x)) x(t) 
/-»oo 

= Xo(0- I dThta{t,T,x)xo(T). (16) 



The kernel of the LHS is bounded and zero for r > f, so it is 
square integrable. Fredholm's theorem therefore guarantees 
the existence of solutions for x. In practice, without exploit- 
ing symmetries or assumptions about the kernel such equa- 
tions are typically solved using numerical techniques [? ? ]. 



Temporal discretization transforms Eq. ( 16 1 into matrix form 



[I-(H„ + H3ct)]:E = [I- 



(17) 



where I is the identity matrix, JE is a measurement record vec- 
tor and the matrices H are the discretized kernels of Eq. ( 16 1. 
The inner product between H and x effects the integration. 
Consequently it is possible to solve for x to finally give 



:E = [I-(H^ + H,,,)]"'[^ 



n] -^o 



(18) 



This expression shows that, for a linear oscillator, it is possi- 
ble to exactly reproduce the measurement record that would 
be obtained with from a system with non-stationary processes 
by causally filtering the measurement record without feed- 
back. This result is valid in both the quantum and classical 
regime and rigorously proves that no force sensing advantage 
is provided by linear feedback onto linear oscillators. This 
precludes enhancements of both sensitivity ifTTl and band- 
width IfTSlI beyond that achievable with estimation alone. 



